Description of the spin structure function gi at arbitrary x and arbitrary Q^ 
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^^ ' The explicit expressions describing the structure function gi at arbitrary x and Q^ are obtained. In 

(w^ \ the first place, they combine the well-known DGLAP expressions for gi with the total resummation 

Cn , of leading logarithms of i, which makes possible to cover the kinematic region of arbitrary x and 

^ ' large Q^ . In order to cover the small-Q^ region the shift Q^ — > Q^ + p? in the large-Q^ expressions 

Q ' for gi is suggested and values of /i are estimated. The expressions obtained do not require singular 

>-^ , factors a::"" in the fits for initial parton densities. 

Q\ • PACS numbers: 12.38.Cy 

(N : 

^: 

JL , The goal of obtaining universal expressions describing the structure function gi at all x and Q^ is an attractive 

ij ' task from both theoretical and phenomenological point of view. Until recently, the only theoretical instrument to 

(~| ] describe gi was the Standard Approach (SA) which involves the DGLAP evolution equationsfl'l and standard fits[2| 

'— '. for the initial parton densities 5q and 5g. The fits are defined from phenomenological considerations at a; ~ 1 and 

Q^ = jj? ^ IGeV^. The DGLAP equations are one-dimensional, they describe the Q^ -evolution only, converting 5q 

and 6g into the evolved distributions Ag and Ag. The DGLAP equations are theoretically grounded in the kinematical 

tH- ' the region A only: 

!>. . A: s>Q2>^2^ ;j<l (1) 

where we have denoted s = 2pq, with p and q being the momenta of the initial hadron and photon respectively. This 
leaves the other kinematical regions uncovered. It is convenient to specify those regions as follows: 
r — ' The small- a; region B: 

B: s:»Q^:»^i'^, x<l (2) 



> 



o 



I. INTRODUCTION 



k>( [ and the small-Q^ regions C and D 
c5 ■ C: 0<Q^<Ai^ a:<l, (3) 



D: 0<Q2<Ai^ x<l. (4) 

As the matter of fact, the SA has been extended from Region A to the small-x Region B, though without any 
theoretical basis. The point is that after converting 5q and 5g into l^q and l^g with the DGLAP evolution equations, 
they should be evolved to the small- a; region as well. The x -evolution is supposed to come from convoluting Ag and 
Ag with the coefficient functions Cdglap- However, in the leading order Cj^Qj^j^p = 1; the NLO corrections account 
for one- or two- loop contributions and neglect higher loops. This is the correct approximation in the region A but 
becomes wrong in the Region B where contributions ~ In (1/a;) are large and should be accounted for to all orders 
in as- Cdglap do no include the total resummation of the leading logarithms of x (LL), so SA requires special fits 
for Sq and 5g. The general structure of such fits (see Refs. [^ ) is as follows: 

Sq = Nx-^ifix) (5) 

where A'' is a normalization constant; a > 0, so a;~° is singular when a; — > and (p{x) is regular in a; at a; —> 0. In 
Ref. pi] we showed that the role of the factor x"" in Eq. (O is to mimic the total resummation of LL performed in 



Refs [3, Q- Similarly to LL, the factor x~°' provides the steep rise to gi at small x and sets the Regge asymptotics for 
gi at x — > 0, with the exponent a being the intercept. The presence of this factor is very important for extrapolating 
DGLAP into the region B: When the factor a;~° is dropped from Eq. ([5]), DGLAP stops to work at a; < 0.05 (see 
Ref. [3] for detail). Accounting for the LL resummation is beyond the DGLAP framework, because LL come from the 
phase space not included in the DGLAP -ordering 



fi^ <kl ^< kl 



< 



<0' 



(6) 



for the ladder partons (k2i _l are the transverse components of the ladder momenta ki). LL can be accounted only 
when the ordering Eq. ^ is lifted and all ki ± obey 



M 



<k-^<{p + qf « (1 - x)2pq w 2pq 



at 



small X. Replacing Eq. © by Eq. ([7]) leads inevitably to the change of the DGLAP parametrization 



a 



DGLAP 



asiQ^ 



(7) 



(8) 



by the alternative parametrization of as given by Eq. (|14p . This parametrization was obtained in Ref. : 6' and was used 
in Refs. [J,l5| in order to find explicit expressions accounting for the LL resummation for gi in the region B. Obviously, 
those expressions require the non-singular fits for the initial parton densities. Let us note that the replacement of 
Eq. ([6]) by Eq. ([7]) brings a more involved /i -dependence of gi. Indeed, Eq. ([6]) makes the contributions of gluon 
ladder rungs be infrared (IR) stable, with /i acting as a IR cut-off for the lowest rung and ki ± playing the role of the 
IR cut-off for the i + 1-rung. In contrast, Eq. ([7]) implies that /i acts as the IR cut-off for every rung. 

The small-Q^ Regions C and D are, obviously, beyond the reach of SA because DGLAP cannot be exploited here. 
Alternatively, in Refs. [3, 0] we obtained expressions for gi in the region C and proved that Region C can be described 
through the shift Q^ — > Q"^ + p? in our large-Q^ formulae. Combining these results with SA obtained in Ref. [3| makes 
it possible to describe gi in Region D. For the sake of simplicity, we present below formulae for gi^ , the non-singlet 
component of gi only. 



II. DESCRIPTION OF gi IN THE REGION B 



The total resummation of the double-logarithms (DL) and single- logarithms of x in the region B was done in 
Refs. 0,0 • In particular, the non-singlet component, g^^ of 51 is 



g^'^ix, Q') - (e^/2) / (1/x)-C^s(^)5q(u;) exp (H^siuj) HQVf^^)) , 



(9) 



with new coefficient functions Cns, 



and anomalous dimensions H 



NS, 



CArs(w) 



^-<*iM 



(10) 



where 



Hns - (1/2) u; - v/w2 _ B{uj) 



B{lu) = (47rCF(l + lu/2)A{lo) + D{uj)) / {2^"^) 



(11) 



(12) 



D{u}) and A{lo) in Eq. (fT^ are expressed in terms of p = ln(l/.T), rj = ln(/i^/AQp^), b = (33 — 2n/)/127r and the 
color factors Cp = 4/3, A^ = 3: 



Diiu) = 



2Cf 



'd^e--ln(^±^) 



p + ri i 

{p + riY + TT^ 77 



(13) 



A{u;) 



dpe-'^P 



if + TT^ Jq [p + ry)2 + 7r2 



(14) 



Hs and Cns account for DL and SL contributions to all orders in as- Eqs. (|14p and p3)) depend on the IR cut-ofF /i 
through variable 77. It is shown in Refs. [J, Q that there exists an Optimal scale for fixing ^: /i « 1 Gev for g^^ and 
/i w 5 GeV for gf . The arguments in favor of existence of the Optimal scale were given in Ref. Q • Eq- & predicts 
that gi exhibits the power behavior in x and Q^ when x —^ 0: 



gr 



{QVx'f^'^', gf 



{Q'h 



1\^sl1 



(15) 



where the non-singlet and singlet intercepts are Ajv5 — 0.42, A5 = 0.86 respectively. However the asymptotic 
expressions P3|) should be used with great care: According to Ref. [3], Eq. pS]) should not be used at a; > 10~^. So, 
Eq. ([9]) should be used instead of Eq. p^ at available small x. Expressions accounting the total resummation of LL 
for the singlet g\ in the region B were obtained in Ref. [a]. They are more complicated than Eq. ([9]) because involve 
two coefficient functions and four anomalous dimensions. 



III. UNIFIED DESCRIPTION OF REGIONS A AND B 



As was suggested in Ref. Q, the natural way to describe g\ in the Regions A and B is to combine the small- a; 
results with the DGLAP expressions for the coefficient functions and anomalous dimensions oi g\. In particular, g^^ 
is again given by Eq. ([9]), however with the new coefficient function Cjqs and new anomalous dimension Uj^s'- 



Cns — Cns ^ 
Hns — Hns 



r<DGLAP A/^ 



NS 



(16) 



nGLAP A TT 

' Ins ~ i^nNS 



where Cats and H^s are defined in Eqs. (jlOIUp . C^g^^^ and ■'Jns^^^ are the DGLAP non-singict coefficient 
function and anomalous dimension. The terms ACns, AHns should be introduced to avoid the double counting. In 
the case when the DGLAP expressions are used in C^'g^"^^ and I'^'s'"^^ with the LO accuracy. 



ACjvs = 1, AHns 



AH 
27r 



(17) 



They are the first terms of expansions of Eqs. pOllip in the series in A{lo). In order to account for the NLO terms for 
^NS^^^ and Jns^^^^ ^^^ ^^^^ terms of the expansions should be included into ACns and AHns- When Eq. (fTB|) 
is substituted into Eq. ^, we arrive at the description oi g^^ covering both Regions A and B. Obviously, the main 
contribution to Cns, Hns at Region A comes from their DGLAP components. On the contrary, the total resumation 
terms dominate at a; ^ 1. When Eq. p6p is used, the initial parton densities should not include singular factors. 



IV. DESCRIPTION OF gi IN THE REGIONS B AND C 



Region C is defined in Eq. ([3]). It involves small Q^, so there are no large contributions In (Q^/ii^) in this region. In 
other words, the DGLAP ordering of Eq. © does not make sense in the region C , which makes impossible exploiting 
DGLAP here. In contrast, Eq. ^ is not sensitive to the value of Q^ and therefore the total resummation of LL does 
make sense in the region C. In Ref. [Tj we suggested that the shift 



Q'+Ai' 



(18) 



would allow for extrapolating our previous results (obtained in Refs. [J, [a] for gi in the region B) into the region C. 
Then in Ref. [g] we proved this suggestion. Therefore, applying Eq. (|18p to g^^ leads to the following expression for 

NO ,. -, . 

5i 



-"-^^ valid in the regions B and C: 



5r(a; + z,Q2) = (e2/2) H ^ (_l_)"c7^5(^)5q(,,) exp (ff^s(^) In ((Q^ + m')^')) , 



(19) 



where z — fi^ /2pq. Obviously, Eq. p^ reproduces Eq. ^ in the region B. Expression for gf looks similarly but more 
complicated, see Refs. [^, Isj for detail. Let us notice that the idea of considering DIS in the small-Q^ region through 
the shift Eq. p^ is not new. It was introduced by Nachtmann in Ref. [lO] and used after that by many authors (see 
6-g- |ll|). being based on different phenomenological considerations. On the contrary, our approach is based on the 
analysis of the Feynman graphs contributing to gi . We also suggest that the following values for /i should be used: 
for the non-singlet component oi gi fi = 1 GeV and fi = 5.5 GeV for the singlet gi. 



V. GENERALIZATION TO THE REGION D 

The generalization of the results of Sect. IV to the Region D can easily be done with replacements 

Cns -^ Cns, Hns -^ Hns (20) 

in Eq. (fT9|) . with Cns, Hns defined in Eq. (|16p . So, we arrive at the final result: the expression for gi which can be 
used in the Regions A,B,C,D universally is 

^(— -) CNs{cj)Sq{u;)exp{HNs{u^)\n{{Q^+fi^)/^,^)). (21) 

We remind that the expressions for the initial parton densities in Eq. (j2ip should not contain singular terms because 
the total resummation of leading logarithms of x is explicitly included into Cns and Hns- 

VI. PREDICTION FOR THE COMPASS EXPERIMENTS 

The COMPASS collaboration now measures the singlet gf at x ^ 10"^ and Q^ < 3 GeV'^, i.e. in the kinematic 
region beyond the reach of DGLAP. However, our formulae for gf^^ and gf obtained in Refs. [3, Q cover this region. 
Although expressions for singlet and non-singlct gi are different, with formulae for the singlet being much more 
complicated, we can explain the essence of our approach, using Eq. (jl9[) as an illustration. According to results of [5:], 
/i « 5 GeV for gf , so in the COMPASS experiment Q^ <^ fp. It means. In {Q^ + fp) can be expanded into series in 
Q^//i^, with the first term independent of Q^: 

g^{x + z, Q\ /.2) = 5f (2, ii^) + Y,{Q'lii^fE^{z) (22) 

fc=i 

where Ek{z) account for the total resummation of LL of z and 

^^—{l/zf[CUu^)Sq{cj) + CUu;)Sgiu;)], (23) 

so that 5q{uj) and 6g{uj) are the initial quark and gluon densities respectively and C|'^ are the singlet coeflficient 
functions. Explicit expressions for C|'^ are given in Refs. [a0|- Therefore, we can makes the following predictions 
easy to be checked by COMPASS: 

A. Prediction 1 

In the whole COMPASS range < Q^ < 3 GeV^, the singlet gi does not depend on x regardless of the value of x. 

B. Prediction 2 

Instead of studying experimental the x-dependence of gf , it would be much more interesting to investigate its 
dependence on 2pq because it makes possible to estimate the ratio Sg/Sq (see Ref. [7,] for detail). 

VII. REMARK ON THE HIGHER TWISTS CONTRIBUTIONS 

In the region B one can expand terms ^ {Q^ + /i^)*^ in Eq. (fT51) into series in (/i^/Q^)" and represent g^^{x + 
z, Q^, iJ?) as follows: 

g^Hx + z,Q\ /.') = g^ix, 0"!^^) + E(A^VQ')'7^fe (24) 

fc=l 



where g^ ^ {x , Q^ / fi^) is given by Eq. ([9]); for explicit expressions for the factors Tk see Ref. Q. The power terms in 
the rhs of Eq. ([M|) look like the power ~ 1/(Q^)'^ -corrections and therefore the Ihs of Eq. ([M)) can be interpreted 
as the total resumniation of such corrections. These corrections are of the perturbative origin and have nothing in 
common with higher twists contributions (= HTW). The latter appear in the conventional analysis of experimental 
date on the Polarized DIS as a discrepancy between the data and the theoretical predictions, with g^ ^ {x , Q'^ / fi^) 
being given by the Standard Approach: 

gNS exp ^ gNSSA _^ jjj.^ (25) 

Confronting Eq. ([^5]) to Eq. ([M]) leads to an obvious conclusion: In order estimate genuine higher twists contributions 
to gf^, one should account, in the first place, for the perturbative power corrections predicted by Eq. ([M)) : otherwise 
the estimates cannot be reliable. It is worth mentioning that we can easily explain the empirical observation made 
in the conventional analysis of experimental data: The power corrections exist for Q^ > 1 GeV^ and disappear when 
Q^ ^ 1 GeV^. Indeed, in Eq. (IM)) ^ = 1 GeV , so the expansion in the rhs of Eq. (IM|) make sense for Q^ > 1 GeV'^ 
only; at smaller Q^ it should be replaced by the expansion of Eq. (fT5|) in ((5^//i^)". 

VIII. CONCLUSION 

The extrapolation of DGLAP from the standard Region A to the small- a; Region B involves necessarily the singular 
fits for the initial parton densities without any theoretical basis. On the contrary, the resummation of the leading 
logarithms of x is the straightforward and most natural way to describe gi at small x. Combining this resummation 
with the DGLAP results leads to the expressions for gi which can be used at large Q^ and arbitrary x (Regions A and 
B), leaving the initial parton densities non-singular. Then, incorporating the shift of Eq. (|18[) into these expressions 
allows us to describe gi in the small-Q^ regions (Regions C and D) and to write down Eq. (PT|) describing gi at the 
Regions A,B,C,D. We have used it for studying the gi singlet at small Q^ which is presently investigated by the 
COMPASS collaboration. It turned out that gi in the COMPASS kinematic region depends on z = fi^/2pq only and 
practically does not depend on x, even at a; <^ 1. Numerical calculations show that the sign of gi is positive at z 
close to 1 and can remain positive or become negative at smaller z, depending on the ratio between Sg and 5q. To 
conclude, let us notice that extrapolating DGLAP into the small-x region, although it could provide a satisfactory 
agreement with experimental data, leads to various wrong statements, or misconceptions. We enlisted the most of 
them in Ref. f!!]. Below we mention one important wrong statements not included in Ref. Q: 

Misconception: The impact of the resummation of leading logarithms of x on the small-x behavior of gi is small. 

This statement appears when the resummation is combined with the DGLAP expressions, similarly to Eq. (J16p . 
and at the same time the fits for the initial parton densities contain singular factors like the one in Eq. ([5]). Such a 
procedure is inconsistent and means actually a double counting of the logarithmic contributions: the first implicitly, 
through the fits, and the second in explicit way. It also affects the small-x asymptotics of gi, leading to the incorrect 
values of the intercepts of gi (see Ref. Q for more detail). 
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